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Anisotropic diffusion processes emerge in various fields such as transport in biological tissue and 
diffusion in liquid crystals. In such systems, the motion is described by a diffusion tensor. For a 
proper characterization of processes with more than one diffusion coefficient an average descrip- 
tion by the mean squared displacement is often not sufficient. Hence, in this paper, we use the 
distribution of diffusivities to study diffusion in a homogeneous anisotropic environment. We de- 
rive analytical expressions of the distribution and relate its properties to an anisotropy measure 
in order to distinguish between isotropic and anisotropic processes. We further discuss the influ- 
ence on the analysis of projected trajectories, which are typically accessible in experiments. For 
the experimentally relevant cases of two- and three-dimensional anisotropic diffusion we derive the 
specific expressions, determine the diffusion tensor, characterize the anisotropy, and demonstrate 
the applicability for simulated trajectories. 
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I. INTRODUCTION 

The random motion of suspended particles in a fluid, 
which is usually referred to as Brownian motion, is an 
old but still fascinating phenomenon. Especially, when 
inhomogeneous [T]-^ or anisotropic media [3HS] are in- 
volved, a lot of questions are still open. From the the- 
oretical point of view, a lot of work has been done [7 
to predict the statistical properties of the trajectories of 
such particles using stochastic methods. On the contrary, 
the development of experiments only recently allows ob- 
taining the paths of individual molecules and particles. 
Especially the observation of two-dimensional trajecto- 
ries using video-microscopic methods, for instance single- 
particle tracking (SPT), is already successfully applied in 
biological systems ^ i9j or to understand the microrheo- 
logical properties of complex Hquids [101 [H]- But also the 
observation of three-dimensional paths becomes feasible 
[T^HH]. The statistical analysis of these trajectories is 
usually accomplished by measuring the mean square dis- 
placement (msd) in order to get the diffusion coefficients 
for the matching theoretical description. However, in the 
anisotropic case the diffusive properties depend on the 
direction of motion and are described by a diffusion ten- 
sor. In such systems, the analysis of msds turned out 
to be not sufficient to determine the anisotropy and ex- 
tract the values of the diffusion coefficients |T5Hl7] . For 
similar reasons, we already introduced the distribution 
of single-particle diffusivities as an advanced method to 
analyze stochastic motion in heterogeneous systems ^18j 
involving more than one diffusion coefficient. 

In the current article, we show the applicability of the 
distribution of diffusivities to analyze trajectories of ho- 



:.tu-chemnitz.de| 



mogeneous anisotropic Brownian motion. We present the 
properties of the distribution as well as their relations to 
established quantities. In order to assess the parame- 
ters of the process, we calculate the characteristic func- 
tion, cumulants and moments of the distribution. For 
the asymptotic decay of the distribution of diffusivities, 
we derive a general expression, which involves the largest 
diffusion coefficient of the system. In conjunction with 
the mean diffusion coefficient of the system, the asymp- 
totic decay enables a distinction between isotropic and 
anisotropic processes. Based on these quantities, we pro- 
vide a measure to characterize the anisotropy of the pro- 
cess from the analysis of SPT data. Since in experiments 
the reconstruction of the complete diffusion tensor is of 
great interest, we extend our concept to tensorial diffu- 
sivities, which offer a simple method to determine the 
entries of the tensor. 

Due to restrictions in SPT experiments the complete 
trajectory is often not accessible [11 [19]. Hence, we in- 
vestigate the influence on the distribution of diffusivities 
and the detection of the anisotropy if only projections of 
the actual trajectory are observed. Even in such cases it 
is possible to estimate bounds of the diffusion coefficients 
from the given projections of the diffusion tensor. Since 
especially two-dimensional and three-dimensional diffu- 
sion processes have a high relevance in experiments we 
apply our considerations to these systems. For homoge- 
neous anisotropic diffusion in two dimensions an analyt- 
ical expression of the distribution of diffusivities exists 
and its moments can be related to the diffusion coeffi- 
cients, which enter the anisotropy measure. Moreover, 
we explain the details to reconstruct the diffusion tensor 
from the tensorial diffusivities as well as from projections 
of the trajectory. Three-dimensional processes are inves- 
tigated analogously although a closed-form expression of 
the distribution of diffusivities does not exist. Addition- 
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ally, we deal with anisotropic processes where one diffu- 
sion coefficient is degenerated corresponding to diffusion 
of uniaxial molecules typical for liquid crystalline systems 

The paper is organized as follows. In Sec.|llj we briefly 
recall the theoretical principles of anisotropic Brownian 
motion based on the diffusion tensor and introduce the 
distribution of single-particle diffusivities, its properties 
and relations to established quantities. To apply our new 
concepts to TV-dimensional homogeneous anisotropic dif- 



fusion processes, we provide in Sec. Ill a general expres- 
sion for the distribution of diffusivities. We demonstrate 
how to distinguish between isotropic and anisotropic pro- 
cesses and explain the reconstruction of the diffusion ten- 
sor. Since in experiments the complete trajectory is often 
not accessible, we investigate in Sec. |IV| how projections 
of the trajectory influence the distribution of diffusivi- 
ties. Finally, in Sec. |Vj we apply our results to specific 
systems of anisotropic diffusion which are typical for ex- 
perimental setups. We substantiate the applicability of 
our findings by analyzing data from simulated anisotropic 
diffusion processes. 



II. DEFINITIONS 

A. Anisotropic diffusion 

An A'^-dimensional anisotropic Brownian motion is 
completely defined by its propagator [21j 



p(x,i|x',t') 



(27r)- 



V[2(t-t')]^dctD 
1 I 



X exp 



-(x-x')^D-i(x-x') 
2 2(t-i') ' ^ ' 



(1) 



where D — O^IDO is the positive definite and symmet- 
ric diffusion tensor, D — diag(i?i,D2, . . . ,Dn) denotes 
its diagonalized form with the diffusion coefficients Di 
belonging to the principal axes, and O is an orthogonal 
tensor which describes the orientation of the principal 
axes relative to the frame of reference. 

For the simulation of such processes an alternative de- 
scription exists, where the propagator corresponds to the 
Langevin equation 



dx 
dt 



(2) 



with %/D = O^VDO and VD = 

dia.g{^/Di,^/D2,- ■ ■ ,VDn)- The vector ^(t) = 
[^i(t), . . . , ^Ar(t)]^ denotes Gaussian white noise in 
N dimensions with = and (<^j(i)0(i')> = 

di,S{t-t')yi,j e{l,2,...,N}. 

Assuming time-translation invariance and considering 
displacements Eq. ([I]) is simplified to p(x' -I- r,T|x') by 



substituting t = t — t' and r = x — x'. This condi- 
tional probability density is averaged by the equilibrium 
distribution po(x') given by the Boltzmann distribution 
to obtain the ensemble-averaged stationary probability 
density 



p(r,r)= yd^x'p(x' + r,T|xXx') 



\ ' exp — r^S^^r 



(3) 



of a displacement r = (ri, . . . , rjv)^ in the time interval r. 
Thus, p(r, r) is an A^-dimensional Gaussian distribution 
with zero mean and covariance tensor S = 2tD. 



B. Distribution of diffusivities 

By observing a trajectory x(t) of an arbitrary stochas- 
tic process in N dimensions individual displacements dur- 
ing a given time lag t can simply be measured for a 
certain particle. Moreover, it is natural to relate each 
displacement to a single-particle diffusivity 



A(r) 



[x(t- 



x(t)] 



2A^t 



(4) 



This simple transformation of displacements to diffusivi- 
ties offers the advantage to compare these quantities for 
different experimental setups and different r. Since for a 
fixed time lag r the single-particle diffusivity is fluctuat- 
ing along a trajectory an important quantity is given by 
the probability density p{D). Therefore, the distribution 
of single-particle diffusivities |TH] is defined as 



p{D,T) = {d [D-Dt{T)]) 



(5) 



where (...) either denotes a time average (...) — 

\miT^oc^/T . . .dt, which is typically accessible by 
SPT, or an ensemble average from different experimental 
methods, such as nuclear magnetic resonance [22]. For 
ergodic systems, as considered here, time average and en- 
semble average coincide. It should be noted that other 
definitions of diffusivity distributions exist in the litera- 
ture [ig. 

For homogeneous systems, i.e., the distribution of dis- 
placements p(r, r) is independent of i, Eq. ([s]) is rewritten 



as 



p(D,t) 



d^rSlD 



2Nt J 



p{r,T) 



(6) 



transforming p(r, t) into the distribution of diffusivities. 

For data from SPT experiments, p(r, t) = {S[r — r(T)]) 
and the distribution of diffusivities is obtained by bin- 
ning diffusivities into a normalized histogram according 
to Eq. ([5]) . The normalization ensures that the histogram 
becomes an estimate of the probability density function. 
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For homogeneous isotropic processes in TV dimensions 
the msd grows hnearly with r, since it obeys the well- 
known Einstein relation (r^(T)) — 2NDcT, where Dc 
is the diffusion coefficient of the process. Due to the 
transformation of displacements to diffusivities by Eq. Q 
the linear dependence on r is removed. Hence, the cor- 
responding distribution of diffusivities becomes station- 
ary and comprises single-particle diffusivities fluctuating 
around Dc- For iV-dimensional homogeneous isotropic 
processes the distribution of diffusivities 



N 
2D, 



exp 



N 

2Dr 



D 



(7) 



is obtained, where r(a;) denotes the gamma function. 
More generally, for normal diffusion in N dimensions the 
Einstein relation holds for large t. In this case, Dc is the 
mean diffusion coefficient of the process. Hence, the dis- 
tribution of diffusivities also converges to the stationary 
distribution given by Eq. ([T]) for homogeneous processes, 
which offers a further advantage of this distribution. 



C. Moments 

The distribution of diffusivities is fully characterized 
by its corresponding moments 



oo 

M„, = {D{tD = JdD D^p{D,t). 



(8) 



It should be noted that the first moment is known as the 
mean diffusion coefficient, which is obtained by a well- 
defined integration. This is in contrast to msd measure- 
ments, where the mean diffusion coefficient is determined 
by a numerical fit to the slope of the msd. By inserting 
Eq. (|6| into Eq. ([s]) the integration of D is accomplished. 
As a result the moments 



M„ 



1 



d^r 



{2NTy 



'p{t,t) 



(9) 



are directly related to the moments of the distribution of 
displacements and, thus, to the moments of the propaga- 
tor p(x, t|x',t'). 



III. PROPERTIES OF THE DISTRIBUTION OF 
DIFFUSIVITIES FOR HOMOGENEOUS 
ANISOTROPIC BROWNIAN MOTION 

A. Distribution of diffusivities 

For homogeneous anisotropic diffusion in N dimen- 
sions, where p(r, r) is a Gaussian distribution with zero 
mean given by Eq. ([s]), the computation of the distri- 
bution of diffusivities, its moments, or its characteris- 
tic function is simplified by reformulating the integral of 



Eq. ([6|. Applying the coordinate transform of r = Qq 

with Q = \/2tO^'\/6 results in the distribution of dif- 
fusivities 



1 ^ \ JV 

i=l / j=l 



x6 \ D 



ofe), (10) 



/2ir 



exp(— Thus, the distribu- 



where P(o,i){qj) 
tion of diffusivities is calculated by integration over in- 
dependent standard normally distributed variables with 
zero mean and unit variance. Since the msd for homo- 
geneous anisotropic diffusion again grows linearly as in 
the homogeneous isotropic case, the r dependency in the 
distribution of diffusivities vanishes. 

By obtaining the distribution of diffusivities, for in- 
stance, from displacements along a single trajectory, in- 
formation about the orientation of the diffusion tensor 
is lost. However, all directions contribute to the distri- 
bution and, thus, it still contains information about the 
diffusion coefficients corresponding to the principal axes, 
i.e., the eigenvalues of D. 



Characteristic function, cumulants and 
moments 



With the transformation Eq. (10), the moments and 



the characteristic function of the distribution of diffusiv- 
ities of anisotropic Brownian motion are calculated. For 
the moments, given by Eq. ([Sl), this yields 



(11) 

So, for instance, the first moment of the distribution of 
diffusivities is given by 



1 ^ 

i=l 



(12) 



which is simply the arithmetic mean of all the diffu- 
sion coefficients Di and coincides with the slope of the 
msd. For higher moments of the distribution of diffusiv- 
ities, it is easier to calculate its characteristic function 
(fc) = (exp(i/ci:i)) = J^AD exp{ikD)p^'^{D) by 



substituting p^ (D) with Eq. (10) and performing the 
Fourier transform to obtain 

^^'(fc) = n / exp (^ifc^ j P(o,i)(9,) 



^ 1 



j=i \/l-ik^ 



(13) 
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From the characteristic function Eq. ( 13 1 the cumu- 
lants of the distribution 



1 9"lnG^^(fc) 



^(m- 1)! 



N 



fc=0 



(14) 

are calculated for m > 0. The moments are recursively 
related to the cumulants by 



M„ 



m— 1 

E 

k=l 



m 



1 



(15) 



It should be noted from Eq. ( 13 1 that the character- 



istic function is a product of independent functions in 
Fourier space. Hence, the distribution of diffusivities 
of an iV-dimensional anisotropic system is determined 
by inverse Fourier transform of the characteristic func- 



tion pg'^iD) = J-i Gg'^(fc) 



tN 



1 ^a/JvC*^) 



where G]^^/j^{k) 



is the Fourier trans- 



form of the one-dimensional distribution of diffusivities 
Po^/NiD) = 1/^/2t:DD,/N exp(-iV£)/(2A)) with dif- 
fusion coefhcient Di/N. Correspondingly, it is obtained 
by convolution of N one-dimensional distributions of dif- 
fusivities 



dAi 



dA 



N 



N 



* P 1 n 



N 



}{D) 



(16) 



which follows directly from Eq. (10 1. Thus, with 



Eqs. (10), (13) and (16), we provide three analogous ex- 



pressions to determine the distribution of diffusivities in 
terms of the eigenvalues of D. Depending on the 
considered system each representation offers its own ad- 
vantages. 



C. Asymptotic decay 

In the following, we present the asymptotic behav- 
ior of the distribution of diffusivities for homogeneous 
anisotropic Brownian motion. We show how it allows 
identifying the anisotropy of the process. 

Considering an M-fold degeneracy of the largest diffu- 



Dm > D 



M+l 



> 



> 



sion coefficient with Di = ■ ■ 
Dn the distribution of diffusivities of the homogeneous 
anisotropic system is obtained from the convolution 

where i^) is the distribution of diffusivities of the 

M-dimensional isotropic system Eq. ([7]) with diffusion co- 
efficient Dc = Di/N, which results from the convolution 
of M identical one-dimensional distributions p^^^^(D). 



For D ^ DiDm+i/ {Di — Dm+i) an approximation is 



introduced and yields the asymptotic behavior of Eq. ( 17 ) 

£) 2 ^ 



pI' 



D- 



N 
2Di 



X exp 



r(f) 




Thus, the general behavior in the logarithmic represen- 
tation is governed by 



logpg^(i?) 



D-i-oo 



N 



(19) 



which corresponds to the asymptotic decay of an A^- 
dimensional, isotropic distribution of diffusivities with 
Doo — max(£'i, 1)2, . . . , -Djv), i.e., an exponential de- 
cay involving the largest diffusion coefficient of the 
anisotropic system. 

In the homogeneous isotropic system describing 
the asymptotic decay is equal to the isotropic diffu- 
sion coefficient Dc, cf. Eq. ([t]), which further coincides 
with the first moment (D). This is in contrast to the 
anisotropic case, where (D) < -Doo- Thus, the devia- 
tion of (D) and Doo offers a measure to characterize the 
anisotropy of the process by 



V 



00 

D) 



(20) 



In cases where both values coincide, i.e., the system is 
isotropic, 77 becomes zero. On the contrary, if one diffu- 
sion coefficient is large compared to the others, {D) — 
Drx,/N resulting in 77 = — 1, which denotes the largest 
possible anisotropy in N dimensions. Thus, rj is a mea- 
sure of the anisotropy, but it is not suitable to compare 
systems of different dimensionality N. It should be noted 
that similar measures with different ranges exist (TSl [H] . 

For experimental data, the decay for large D is ob- 
tained from a fit to f{D) — cexp(— Afit-D) to calculate 
Doo = -^/(2Afit). The actual dimensionality M of pro- 
cesses observed in iV > AI dimensions is estimated with 
M = 2 (D) Afit and, thus, the anisotropy is 77 = N/M — 1. 
For example, if an observed iV-dimensional motion yields 
the largest anisotropy value of 77 — iV — 1, the process is 
effectively a one-dimensional motion. 



D. Reconstruction of D 

For experiments it is of great interest to reconstruct 
the diffusion tensor D from the measurements. If com- 
plete information about the trajectories is available, the 
diffusion tensor of the homogeneous anisotropic process 
can be estimated via the displacements. By defining ten- 
sorial diffusivities analogously to Eq. ^ 

[x, (t + r) - X, (t)] [xj {t + T)-Xj (t)] 



2t 



(21) 
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where Xi{t) denotes the i-th component of the N- 
dimensional trajectory x(t), the hnear r dependence of 
the mixed displacements is removed. These tensorial dif- 
fusivities are simply averaged 



and the eigenvalues of Df^ } are 
<ra) > <ri) ^ • ■ ■ ^ Va e {1, • . ■,N}, 



D, 



(22) 



providing an estimator for the corresponding elements 
of D. Here, (...) either denotes a time average or an 
ensemble average depending on the available data. 



IV. PROJECTION TO AN A/-DIMENSIONAL 
SUBSPACE 

Due to experimental restrictions the complete trajec- 
tory is often not accessible but its projection on an M- 
dimensional subspace can be measured. Such processes 
are commonly known as observed diffusion [251 121] • 

The projection of the distribution of displacements 
Eq. ([3| on the considered subspace is the marginal prob- 
ability density 



p(r^,r)= Idr, 



dra„_^, p(r,r), (23) 



where Tq. ,z = 1, . . . , {N — M) denotes (TV — M) arbi- 
trarily chosen directions which are integrated out. The 
vector a = {a.i . . . un-m) contains the indices ai de- 
scribing which elements of r are omitted. Alternatively, 
the projected distribution of displacements is computed 
by the M-dimensional inverse Fourier transform of the 
characteristic function of p{r, r) where the components, 
= 0, Vi € {1, . . . ,N — M}, which correspond to the 
chosen directions, are discarded. Hence, the distribution 
of displacements of the subspace is 



P(r^^r) 



1 



(27r 



exp 



-i(k 



A/- 



Tj.Af 



G(k, 



(24) 

with the characteristic function of the projected propaga- 
tor G(k*^) = exp[-(k^)Ts^^k*^] . The vector k^^ is an 
A/-dimcnsional sub-vector of the complete fc-space and 
denotes a principal M x M submatrix of S obtained 
by deletion of rows and columns with same respective 
index a^. 

The distribution of diffusivities of such a projected dif- 
fusion process is calculated analogously to Eq. ^ by 
integrating over r^^. Since S = 2rD is a symmetric, 
positive definite matrix for r > 0, all principal submatri- 
ces S*^ are symmetric, positive definite matrices as well. 
Hence, they can be diagonalized and, thus, the projected 
distribution of diffusivities has the Af-dimensional form 
of the generic expression Eq. (10), or (16) 



ever, it depends on the eigenvalues 

of the projected diffusion tensor D*^ 
eigenvalues of D are identified as 



-lA/ 



1, 



How- 
. . , M 



Di> D2> ■ ■ ■ > D 



N 



T,'^/{2t). If the 



(25) 



the well-known interlacing inequalities [17] require 



(26) 



Dk > Dl^^^^ > Dk+u VA; e {1, . . . , iV - 1} (27) 

for all a G {1. ... , N}. This expression is applied recur- 
sively (N-M) times to obtain a relation for the eigenvalues 
of the principal M x M submatrix [28] 

Dk > Dt'^ > Dk+N-M, Vfc e {1, . . . , M} (28) 

for arbitrary a. By implication, if at least two eigenval- 
ues of the submatrix D*^ differ, i.e., the projected process 



is anisotropic, Eq. (27) states recursively that the com- 



plete process is anisotropic as well. Thus, the distribution 
of diffusivities of the projected iV-dimensional anisotropic 
Brownian motion may already indicate the anisotropy of 
the complete process as well as the magnitude of one of 
the involved diffusion coefhcients. However, a single pro- 
jection is not sufficient to obtain the underlying diffusion 
coefficients. 

Nevertheless, it is possible to estimate the bounds of 
the diffusion coefficients. The lower bound of the eigen- 
values is given by zero, due to the positive semidefinite- 
ness of D. An upper bound for the largest eigenvalue 
can be found if enough projections or submatrices are 
available to cover the whole trace of D. By use of the 
relation between the trace of an x iV matrix A and its 
eigenvalues Ai, tr A = subtotals of the trace of D 

are given by the sum of the eigenvalues of the respective 
submatrices. If the non-overlapping orthogonal projec- 
tions of D defined by a compose a partition of the set 
{1, . . . , N}, the trace of the tensor is given by 



N 



(29) 



with IJa = {1, . . . , N}, where the partition elements a 
do not necessarily have identical dimensionality. 

For example, if one measures the eigenvalues of two 
non-overlapping projections of a 3 x 3 diffusion tensor D, 
the trace of D is given by 



trD = trDj^ 3) + trD(2) 
= ^1,(1 3) + ^h{2) + ^ 



2 

2,(2)- 



(30) 



Thus, the eigenvalue inequalities for that example using 
the relations above are given by 

trD >Di > max(i:ii ^ 3 ( _^)) > D2 

D2 >min(i?i\(, 3),i^2,(2)) > ^3 > 0, (31) 

which allows a rough estimation of the diffusion coeffi- 
cients from the given projections. 
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V. SPECIFIC SYSTEMS 

A. Two-dimensional systems 

The distribution of difFusivities of a two-dimensional 
homogeneous anisotropic system is calculated, for in- 
stance, via Eq. (16), resulting in 



oo oo 

= y dAi y dAs 5[D- (Ai + A2) 




exp 



2 I Di 



1 

£>2 



D 









[I 




D2) 



D 

(32) 



where Iq{x) denotes a modified Bessel function of the 
first kind. The first two moments of this distribution, as 
given by Eq. (|8l, yield 



{D)=Ah = ^pi+Z?2) 



(33) 




Figure 1. The distribution of diffusivities (histogram) from a 
simulated trajectory of a homogeneous anisotropic diffusion 
process in two dimensions with diffusion tensor D given by 
Eq. ( 39 1 agrees well with the analytic distribution of diffu- 
sivities (solid line) from Eq. ( |32[ ) with Di = 5 and D2 = 1 
denoting the eigenvalues of D. Additionally, the asymptotic 
function Eq. (36 1 (dotted line, D^o ~ 5) agrees reasonably for 



and 



large D. Furthermore, a distribution of diffusivities (dashed 
line) of two-dimensional isotropic diffusion with the same av- 
erage diffusion coefficient Dc = (D) — {Di + D2)/2 = 3 is 
shown for comparison. 



(D"^) = Ah = -.{iDl + 2DiD2 + iDl) 



(34) 



Hence, the mean diffusion coefficient coincides with the 
arithmetic mean of the diffusion coefficients belonging to 
the two directions of the anisotropic system as expected 
from Eq. (12). Solving the simultaneous Eqs. (33) and 



( 34 ) yields the expression 



Di,2 = Ml ± VM2 - 2M^^ 



(35) 



to obtain the diffusion coefficients Di and D2 from the 
moments. 



The asymptotic behavior of Eq. (32) for large D is 
given by Eq. ( 18 ) and yields 



exp 



D 



^\D^-D2\ttD 



(36) 



max(Z?i, 1)2). Thus, the behavior in the 



with i'oc 

logarithmic representation is governed by 



D- 



D 



(37) 



which corresponds to the decay of the distribution of 
difFusivities in two-dimensional homogeneous isotropic 
systems with 13 00, i-e., an exponential decay with the 
largest diffusion coefficient of the anisotropic system. Ac- 
cordingly, the smallest diffusion coefficient is given by 
2 (D) — Dao- From the asymptotic decay and the mean 



diffusion coefficient the anisotropy of the system is char- 



acterized by Eq. ( 20 ) and corresponds to the ratio 



ri = 



\Di-D2\ 
D1+D2 



y/M2 - 2Mf 
Ml 



(38) 



which is also related to the moments. 

To substantiate our analytical expressions by results 
from simulation, a random walk was performed by nu- 
merical integration of the Langevin equation Eq. ^ in 
two dimensions using the diffusion tensor 



D 



4 

V3 2 



(39) 



with the eigenvalues Di = 5 and D2 — I. The obtained 
trajectory of two-dimensional homogeneous anisotropic 
diffusion consisted of 10^ displacements and its distri- 
bution of diffusivities is depicted in Fig. [T] The agree- 
ment of the normalized histogram from simulated data 
with the analytic distribution Eq. (32) is obvious. Devia- 



tions between simulation and the analytic curve for large 
D are due to insufficient statistics from the finite num- 
ber of displacements. Moreover, Fig. [T] shows the mono- 
exponential behavior corresponding to isotropic diffusion 
in two dimensions for comparison. Although the average 
diffusion coefficients of both processes coincide, the de- 
cays of the distributions differ. The reason is the asymp- 
totic behavior given by Eq. ( 36 ) in the anisotropic case 



which exponentially decays with the largest eigenvalue 
for large D as depicted in the figure. In contrast, for 
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the isotropic system the asymptotic decay corresponds 
to the mean diffusion coefficient resulting in the quan- 
titative difference. Furthermore, the distributions are 
quahtatively different for small D. A characteristic dif- 
ference between the systems is the convex shape in the 
logarithmic representation of the anisotropic distribution 
of diffusivities. This intuitively results from the two dif- 
ferent exponential decays related to the distinct diffusion 
coefficients Di and D2. In a more rigorous way, since 
3752 logp^(-D) > 0, the anisotropic distribution of diffu- 
sivities is a superconvex function [29| . 

For experimental data, it is easy to calculate the first 
two moments Mi and M2 by averaging the short-time dif- 
fusivities of Eq. Q and their squares, respectively. The 
averaging is accomplished either along a single trajectory 
or from an ensemble of trajectories avoiding any numer- 
ical fit. The first two moments are sufficient to calculate 



Di and L»2 by Eq. (|35|). 

For the sample trajectory used in Fig. [T] the first 
two moments are determined to be Mi = 2.987 and 
M2 = 21.72. According to Eq. (35), the underlying dif- 
fusion coefficients yield Di — 4.956 and D2 = 1.018. 
These values agree well with the eigenvalues of the ten- 
sor Eq. ( 39 1 , which was used as input parameter of the 
simulation. The calculation of 77 = 2/3 indicates a con- 
siderable anisotropy of the process. 



Limiting cases 



In the case of identical diffusion coefficients for both di- 



rections the anisotropy vanishes as discussed for Eq. ( 38 1 



The resulting isotropic diffusion process is character- 
ized by a single diffusion coefficient Dc = Di = D2- 



Hence, Eq. (32) simplifies to the well-known distribution 
of single-particle diffusivities of two-dimensional isotropic 
diffusion m 



exp 



(-*) 



(40) 



given by an exponential function. 

If, on the contrary, the anisotropy is large, diffusion in 
one direction will be suppressed. Without loss of gener- 
ality, this is accomplished by tending one of the diffusion 
coefficients to zero. Thus, by taking the limit of vanishing 
D2, the distribution of diffusivities Eq. (32) is simplified 
to 



exp 



D 

Di 



(41) 



which has the structure of the distribution of diffusiv- 
ities of one-dimensional diffusion |18) . Since diffusion 
into the perpendicular direction is prohibited, as ex- 
pected, it qualitatively leads to the observation of a one- 
dimensional process. This can be identified by the char- 
acteristic factor D~^/^ due to which the distribution of 



diffusivities diverges for small D. Applying Eq. ([s]) the 
first moment of Eq. (41), i.e., the mean diffusion coeffi- 
cient, yields (D) = DYj2. The factor of 1/2 results from 
the single-particle diffusivities Eq. ^ with N — 2 assum- 
ing that a two-dimensional process is observed. However, 
due to the suppression of one direction this assumption is 
no longer valid and = 1 should have been used instead. 
This conclusion is also obtained from the anisotropy value 
T] ~ 1, which is equal to its maximum value for two- 
dimensional anisotropic processes since effectively one- 
dimensional motion is observed. 

The alternative way to attain a large anisotropy would 
be to prefer diffusion in one direction by increasing the 
corresponding diffusion coefficient. Since there is always 
a contribution from the non-vanishing perpendicular di- 
rection, the distribution converges to {DiD2)'^^^ for 
small D. In the limit of an infinite diffusion coefficient, 
i.e., Di ^ D2, a rescaling of D with (D) keeps the first 
moment constant and finally results in a behavior iden- 
tical to Eq. (fill). 



2. Reconstruction o/D 

In addition to the eigenvalues, it is sometimes of inter- 
est to determine the orientation of the principal axes of 
the system relative to the given frame of reference. This 
is achieved by the reconstruction of the diffusion tensor 



D 



Dii D12 
D12 D22 



(42) 



where the off-diagonal elements are labeled identically 
due to symmetry reasons. The reconstruction is accom- 
plished in two ways either by considering the complete 
two-dimensional trajectory or by using one-dimensional 
projections of the trajectory. 

In the first approach the tensorial diffusivities of 
Eq. ( 21 ) are used to obtain the tensor entries of D. In ac- 



cordance with Eq. ( 22 ) the tensor elements are estimated 



by averaging the tensorial diffusivities along a trajectory 
or over an ensemble. Moreover, the eigenvalues of tensor 
D are expressed by its entries 

1 



Di^2 = 2 (Dii + D22 ± ^{Dii - D22Y + ^Dl^j (43) 

and correspond to the diffusion coefficients of the system. 

For the sample trajectory used in Fig. [T] the measured 
values Dij yield the diffusion tensor 



D 



3.983 1.719 
1.719 1.990 



(44) 



which agrees reasonably with the input parameters of the 
simulation. The eigenvalues from this measured tensor 
Di = 4.973 and D2 = 1.000 show a good agreement with 
the eigenvalues of the input tensor Di — 5 and D2 = 1. 

The second approach determines the tensor D exclu- 
sively from one-dimensional projections of the trajectory. 
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In order to obtain results, at least three different projec- 
tions are necessary. For simplicity, it is preferable to use 
projections along two perpendicular axes, which define 
the frame of reference for D. Further, a projection onto 
an axis is required which is rotated about an angle 6 rel- 
atively to the frame of reference. In such a setup, the 
first moments of the distribution of diffusivities related 
to the first two projections are identical to the averaged 
tensorial diffusivities (DI^{t)) and (Z?p(r)). Thus, they 
yield the two diagonal elements of D. The first moment 
of the third projection measures the leading diagonal cl- 
ement Df 1 of the rotated tensor D" = R(0)'^DR(6') with 

rotation tensor R(0) = I ^^"^^ 
^ ' y^smt' cost' 

value is sufficient to obtain the off-diagonal element of D 

from 



This additional 



-Dl2 = 



D^i - cos2 B - D22 sin^ 
sin(26') 



(45) 



For the calculation, any projection of the trajectory onto 
an arbitrary one-dimensional axis, i.e., any 0, can be used 
except directions perpendicular or parallel to axes of the 
frame of reference, i.e., angles 9 which are multiples of 
7r/2. It should be noted that the reconstruction from 
the distribution of diffusivities of projected trajectories is 
possible although the definition of the diffusivities omit 
any directional information. 

In the example with Dn — 3.983, D22 = 1.990 and a 
measured U^^^^ — 2.983, the off-diagonal element yields 
D12 = 1.719, which is in good agreement with the value 
« 1.732 appearing as input parameter of the simula- 
tion. 

In conclusion, it depends on the constraints of the ex- 
periment which of both approaches is more practicable. 
In either way the complete diffusion tensor D is recon- 
structed reasonably. 



B. Three-dimensional systems 



Analogous to the two-dimensional case, it is possible to 
calculate the distribution of diffusivities either by inverse 
Fourier transform of the general characteristic function 
Eq. ( 13 1 or by the convolution Eq. ( 16 ). In both cases the 



analytical integration cannot be performed completely. 
However, the integration can be accomplished numeri- 
cally. By integrating two variables Eq. ( 16 ) is reduced 



to 



00 00 



p'^iD) = / dAi^ dA2 j dAg 5 [D - (Ai + A2 + A3) 





3 ' 



dAi 



3/2 



1 







X exp 

X /o 



2 J ^ttDiD2D3Ai 
1 1 



3 

4 [\D2 D3 
3/1 1 

1\D3'D2 



{D-A,) 



2Ai 



(46) 



For further simplification, a series expansion of the mod- 
ified Bessel function /o(a;) can be applied, which allows 
performing the last integration. However, this only re- 
sults in a converging sum, which cannot be simplified any 
further. 

By using the general expression of the cumulants 
Eq. ( 14 1 and the relation between cumulants and mo- 



ments Eq. (15), the first three moments of the distri- 



bution of diffusivities of three-dimensional homogeneous 
anisotropic diffusion processes are 



1 



Ml = -iDi+D2 



(47) 



M2 
and 



[{Di 



D2 + D3f + 2{Dl + Dl + Di)], (48) 



Af3 = - [5Dl + 'SDl{D2 + D3) 



' 2D2D3 



-{D2 + D^){bDl - 2D2D^ + bDl)] . (49) 



These expressions are similar to Eqs. (33) to (34) and re- 



late the moments of the distribution of diffusivities to dif- 
fusion coefficients Di to D3 of the anisotropic process. By 
solving the simultaneous Eqs. (47) to (49) the underly- 



ing diffusion coefficients are determined by the measured 
moments of the distribution. The solution comprises six 
triplets {Di to D3), which are permutations of the three 
diffusion coefficients. Due to the cubic contributions in 



Eq. (49) the expressions are too large to show here but 



can be easily obtained. 

The asymptotic behavior of Eq. (46) for large D is 
given by Eq. (18), which assumes Di > D2 > D^, and 
results in 



D- 



^^exp(^^) 
^27r{D, - D2){D, - D3)D' 



(50) 



Thus, the behavior in the logarithmic representation is 
determined by 



D^c 



3D 

2DZ. 



(51) 
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Figure 2. Distribution of difFusivities (histogram) from a sim- 
ulated trajectory of a homogeneous anisotropic diffusion pro- 
cess in three dimensions with diffusion tensor D given by 



Eq. ( 53 1 agrees well with the distribution of diffusivities (solid 
line) obtained from numerical integration of Eq. ( |46[ ), using 
the eigenvalues Di — 5, D2 = 3 and D3 = 1 of tensor D. 
For comparison, the distribution of diffusivities (dashed line) 
of an isotropic diffusion process in three dimensions, given by 
Eq. (54 1, is shown with the same average diffusion coefficient 
Dc = {D) = {Di + D2 + £'3)73 = 3 as in the anisotropic pro- 
cess. In the inset, the asymptotic function Eq. (50 1 (dotted 
line, Di = 5) agrees reasonably for large D. 



case the crossover in the logarithmic representation ex- 
hibits a concave curvature while in the anisotropic case it 
is convex after the peak. This is expected due to the dif- 
ferent exponential decays related to the distinct diffusion 
coefficients Di to analogously to two-dimensional sys- 
tems. However, a rigorous proof of ^^logp^(r') > 
for D > Dpeak is still missing. In addition, the asymp- 
totic behavior given by Eq. ( 50 ) is depicted and provides 



a reasonable approximation for large D. 

The eigenvalues of D for experimental data arc easily 
determined by measuring the leading moments of the dif- 
fusivities. For the sample trajectory used in Fig. [2] the 
first three moments result in Mi = 2.995, M2 = 16.68 
and M 3 — 140.2. By solving the simultaneous Eqs. (47) 



to ( 49 1 , the underlying diffusion coefficients yield Di 



4.884, D2 = 3.153 and D3 = 0.948. These values agree 
well with the eigenvalues of the tensor Eq. (53), which 



was used as input parameter of the simulation. The cal- 
culation of 77 = 2/3 indicates a considerable anisotropy 
of the process. 



1. Limiting cases 

If the diffusion coefficients of all three directions coin- 
cide with Dc — Di = D2 = -D3, the distribution of diffu- 
sivities for the three-dimensional isotropic system |l8j 



which corresponds to the asymptotic decay of a three- 
dimensional isotropic distribution of diffusivities with 
Z?oo = mBx(Di, D2, D3). The anisotropy parameter 



Eq. ( 20 ) in the three-dimensional case corresponds to 



V = 



(Di - D2) + {Di - D3) 
D1+D2 + D3 



(52) 



which considers the differences of the individual diffusion 
coefficients to characterize the anisotropy. It is obvious 
that the largest anisotropy yields j] = 2. 

In order to substantiate our results by simulated 
data, the simulation of a three-dimensional homogeneous 
anisotropic random walk was performed using the diffu- 
sion tensor 



D 



V 



4 _i/3 _i\ 
^ 22 

\/3 13 3^3 
2 4 4 

_ 1 3\/3 

2 4 4 



(53) 



The obtained trajectory consists of 10^ displacements 
and its distribution of diffusivities is depicted in Fig. [2] 
The distribution of diffusivities from the simulated tra- 
jectory shows a good agreement with the curve obtained 
from numerical integration of Eq. (46). The deviations 



for larger values of D result from the finite simulation, 
i.e., its insufficient statistics. Furthermore, Fig. [2] shows 
the distribution of an isotropic system where a qualitative 
distinction at the crossover from the maximum peak to 
the exponential decay becomes apparent. In the isotropic 



(54) 



will be obtained from Eq. ( 46 ) in agreement with Eq. ([7| . 

If exactly two diffusion coefficients coincide, one usu- 
ally refers to diffusion processes of uniaxial molecules 
|20j . In this case, the general distribution of diffusivities 
of three-dimensional homogeneous anisotropic diffusion 
Eq. ( 46 ) simplifies to 



3 exp I 
2 



3D 
2_D(2) 



erf 



(\/2(_D(i) 



D 



^£,(2) (£,(2) _£)(!)) 



(55) 

where D'-^'> and D'-^^ are the eigenvalues of D with mul- 
tiplicity one and two, respectively. In general, a dis- 
tinction between the oblate case (D^^^ > D'-^\ disc) 
and the prolate case {D^^^^ < D^^\ rod) is drawn for 
uniaxial molecules. In the prolate case both square 



roots in Eq. (55) yield complex numbers. However, with 
erf {^—x) I \f—y = erfi( y^) / ^ for x, y > and x^y € R, 
Eq. ( |55[ ) remains a real valued function. Hence, a distinc- 
tion between the two cases of the diffusion coefficients is 
not required for the distribution of diffusivities. 

In the uniaxial case the first three moments simplify 

to 



1 



Ml = -(£»(!) +2i:>(2)) 



M2 = \ (3i?(i)' + 47^(i)l?(2) + 8i^(2)2)j 



(56) 



(57) 
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and 



9 V 

Thus, the eigenvalues of D are calculated by 



(58) 



= Afi T V3Af2 - 5M2 



and 



(59) 



(60) 



where the sign in the equations depends on the constraint 
of positive diffusion coefficients. None of the eigenvalues 
wiU become complex since with Eqs. (56 1 and (57) the 
expression under the square root 3M2 — 5Af2 = | (S'-^' — 
£)(2))2 > is always positive and, hence, M2 > fMf. ft 
should be noted that for f Mf < M2 < 2Mf both signs 



in Eqs. (59) and (60) yield positive diffusion coefficients. 
In this case, the third moment has to be exploited in 
order to decide the correct pair of diffusion coefhcients 
by comparing Eq. (58 1 with the measured value. Hence, 



there exist distributions of diffusivities with identical mo- 
ments Ml and M2, which result from different diffusion 
coefficients. In this case, the distinct M3 determines the 
corresponding diffusion coefficients of the system. In the 
limit M2 ^M^, D^^) and £'(2) approach each other. 
In this particular case, the decision for the correct pair 
cannot be made accurately since both pairs yield approx- 
imately the same M3 from Eq. (58). However, this limit 



corresponds to the isotropic system and, hence, the single 
diffusion coefficient is directly given by the first moment 
of the distribution. 

Fig. |3] depicts examples of such distributions for the 
general anisotropic, the prolate, and the oblate case. The 
differences can be identified qualitatively. In the gen- 
eral and in the prolate case, the decay after the max- 
imum peak has a convex curvature in the logarithmic 
representation, whereas in the oblate case it decays in a 
purely concave manner. This qualitative change is ob- 
tained from logp^^i(il'), where for I?(2) > £)(i) this 
second derivative is always negative, i.e., p'^'^^{D) is a su- 
perconcave function. In the opposite case a zero-crossing 
is observed which results in a change from concave to 
convex curvature. However, a rigorous proof for the lat- 
ter is still missing. In all cases, the exponential decay for 
large D is determined by the largest diffusion coefficient 
as given by Eq. ( 50 1 . However, since the first decay after 



the peak is dominated by the smallest diffusion coeffi- 
cient, the curve is shifted to the left for the prolate case 
in contrast to the oblate case when D2 is changed from 
D3 to Di. As expected from the first moment, the gen- 
eral case lies in between. A better distinction between the 
different cases is achieved quantitatively by determining 
the moments and calculating the diffusion coefficients. 

In Fig. |4] the distribution of diffusivities for different 
ratios 




D 



Figure 3. Distribution of diffusivities (lines with open sym- 
bols) of different homogeneous anisotropic diffusion processes 
in three dimensions. A qualitative distinction between the 
oblate case (O; D^^' = l,D^^^ = 5), the prolate case (O; 
-D'^-* = 5,D'^^ — 1), and a general anisotropic case (O; 
Di = 5, D2 = 3, -D3 — 1) is possibie since the decay after the 
maximum peak shows a concave curvature in the first case 
and a convex curvature in the latter cases. Furthermore, each 
anisotropic case obeys the same asymptotic decay given by 
the largest diffusion coefficient. For comparison the isotropic 
cases (lines with filled symbols) with Dc = 1 (■) and Dc — 5 
(a) are given, which always have a concave shape and, thus, 
are qualitatively indistinguishable from the oblate case. 



is shown, ranging from oblate cases (r < 1) to prolate 
cases. It can be seen that in the limit D^^^^ — > and, thus, 
r — )• 0, the distribution converges to the two-dimensional 
isotropic case with Dc — 2/3Z?^2) Por r — 1, the distri- 
bution converges to the three-dimensional isotropic case. 
In the prolate cases the distribution separates signifi- 
cantly from the three-dimensional isotropic case for in- 
creasing r. For further increasing ratios (r — )■ 00) the 
distribution converges to the one-dimensional isotropic 
case with Dc — 1/3D^^\ In contrast, the oblate cases 
converge rapidly to the two-dimensional isotropic case 
for decreasing r. A qualitative distinction may only be 
possible for small D, where the distribution still deviates 
from the mono-exponential behavior of the isotropic sys- 
tem. Howev er, quantitatively the anisotropy is character- 



ized by Eq. (|20| , which results in 77 = and rj = ^^2+r^ 



D^'^/D 



(2) 



(61) 



for oblate and prolate cases, respectively. Thus, in the 
oblate case the largest possible anisotropy emerges at 
small r, which yields 77 = 1/2 and clearly indicates the 
anisotropy. In the prolate case, the largest anisotropy 
will be obtained, if only one direction is preferred. Then, 
the anisotropy parameter 77 = 2 is maximal, which corre- 
sponds to one-dimensional motion in a three-dimensional 
system. 
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10' 



r = 0.1 

r = 0.5 — 0- 

r = 2.0 --0- 

r= 10.0 --E3-- 




of the diffusion coefficients. As explained in Sec. IV any 



Figure 4. Distribution of diffusivities ( line s with open sym- 
bols) for different ratios r given by Eq. (61 1 and fixed D^^'' = 
1. The crossover from oblate cases (r < 1, solid lines) to 
prolate cases (r > 1, dashed lines) shows a broadening of the 
peak for increasing ratios. Again, the behavior after the peak 
changes from concave to convex, respectively. For compari- 
son, the distribution of diffusivities of the limiting isotropic 
cases are depicted for two-dimensional (■, Dc — 2/3) and 
three-dimensional processes (A, Dc — 1). The distinction of 
prolate cases from the isotropic limits is simpler than for the 
oblate cases. 



2. Reconstruction of D 

As discussed for two-dimensional processes, the dif- 
fusion tensor will be easily obtained by measuring the 
averaged tensorial diffusivities according to Eq. (22) if 



the complete three-dimensional trajectory of the homo- 
geneous anisotropic process is available. For the sample 
trajectory used in Fig. [2] the measured values Dij yield 
the diffusion tensor 



D 



3.994 
-0.862 
-0.492 



-0.862 
3.233 
1.296 



-0.492> 
1.296 
1.758 



(62) 



which agrees reasonably with the input parameters of the 
simulation Eq. (53). Further, the eigenvalues from this 
measured tensor Di = 4.983, D2 = 2.998 and D3 = 1.004 
show a good agreement with the eigenvalues of the input 
tensor Di = 5, I?2 = 3 and = 1. 

However, if only a projection of the complete trajectory 
is available, e.g. from SPT, only the properties of the re- 
spective submatrix of D can be measured. For instance, 
if the two-dimensional projection onto the x-y-plane of 
the sample trajectory is available, the first two moments 
of the distribution of diffusivities are determined to be 
Mi_^ = 3.613 and M2,z = 26.95. Using Eq. ([sH]), the 
eigenvalues of the principal submatrix are computed 
to he Dj^ = 4.531 and I)| ^ = 2.695. Hence, the eigen- 



value inequalities of Eq. ( 27 ) provide the estimate 



> Di 



4.531 >D2>Dl^ 



2.695 > > 0. 



further observed projection improves the estimates of the 
eigenvalues of D. An additional projection onto the x-z- 
plane, for instance, yields the moments Mi^y — 2.876 and 
^h,y — 18.00 resulting in the eigenvalues D\ y — 4.083 
and D^^y = 1.668. Since with two orthogonal two- 
dimensional projections of the three-dimensional process 
the complete trace of D is covered, an upper bound 
for the largest eigenvalue is found to be Di < trD < 



n2 



Dl y Dl y = 12.977. Hence, the eigen- 



value inequalities yield 

12.977 >Di> max{Dl^,Dly) = 4.531 
min{Dl^,Dl y) = 4.083 > D2 > max(i)^^^, L>|j^) = 2.695 
min{Dl^,Dly) = 1.668 > i^s > 0. (64) 

If additionally the projection onto the y-z-plane is avail- 
able the eigenvalues of D arc estimated more precisely 
similar to the previous steps. To improve the upper 
bound of Di , the trace of D is calculated from all these 
eigenvalues by trD = i(^2^ ^ f)2^^ ^ f)2^^ _^ ^2^^ ^ 

Df ^ + DI z), where the prefactor arises from the over- 
lapping traces of the submatrices. 

In case of the availability of all orthogonal two- 
dimensional projections of the process the tensorial dif- 
fusivities offer an advanced approach to determine the 
diffusion tensor. Since their first moments yield the en- 
tries of the principal submatrices D2, T>1 and T>1 the 
underlying diffusion tensor D is completely defined. 

To summarize, the experimental setup influences the 
available data and affects how many parameters of 
the underlying process can be restored. A single 
two-dimensional projection may already hint at the 
anisotropy of the process. However, it is not sufficient to 
give an upper bound for the largest eigenvalue. An ad- 
ditional orthogonal two-dimensional projection or even a 
one-dimensional projection in the missing direction de- 
termines this upper bound and narrows the ranges of the 
eigenvalues. For a reconstruction of the complete tensor 
either the complete trajectory or three orthogonal two- 
dimensional projections of the process are necessary. 



(63) 



VI. CONCLUSIONS 

To investigate A'^-dimensional homogeneous 
anisotropic Brownian motion we applied the distri- 
bution of diffusivities obtained from single-particle 
tracking data. From the general expression of the 
distribution of diffusivities we derived the relation of 
its moments and cumulants to the eigenvalues of the 
diffusion tensor D. We showed that the asymptotic 
behavior of the distribution of diffusivities deviates 
from the isotropic case which allowed establishing an 
anisotropy measure based on the distribution. Since, 
due to experimental restrictions, often only projections 
of the trajectories are observed we further discussed the 
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consequences and provided an estimate for the bounds 
of the involved diffusion coefficients. 

After our general considerations, we applied the results 
to specific systems with a high relevance in experiments. 
In particular, we investigated two-dimensional and three- 
dimensional systems as well as uniaxial molecules in 
three dimensions. In a two-dimensional homogeneous 
anisotropic system, the distribution of diffusivities com- 
prises a modified Bessel function and allows a qualitative 
distinction from the mono-exponential decay observed in 
isotropic systems. Moreover, the first two moments of 
the distribution are sufficient to calculate the diffusion 
coefficients corresponding to the principal axes. Even 
the orientation of the principal axes and, thus, the com- 
plete diffusion tensor D can be determined by using ten- 
sorial diffusivities or via three one-dimensional projec- 
tions of the trajectory. For three-dimensional processes 
the general expression of the distribution of diffusivities 
is more elaborated and the last integration has to be 
evaluated numerically. However, we expressed the first 
three moments in terms of the diffusion coefficients be- 
longing to the principal axes. Conversely, these expres- 
sions offer a method to calculate the diffusion coefficients 
from the moments measured in experiments, where other 
analysis fails. It is further shown that the isotropic and 
anisotropic systems differ in the logarithmic representa- 
tion of the distribution of diffusivities, i.e., the asymp- 
totic decay rate is proportional to the inverse slope of the 
msd and to the inverse of the largest diffusion coefficient, 
respectively. Thus, the distribution of diffusivities for 
anisotropic diffusion asymptotically decays slower than 
for isotropic diffusion with the same average diffusion co- 
efficient. The deviation between the asymptotic decay 
and the decay belonging to the first moment provides a 
suitable measure for the anisotropy of the process. For 
uniaxial molecules diffusing in three dimensions the third 
integration is accomplished and the resulting distribution 
of diffusivities involves an error function. In this case, the 
diffusion coefficients of the principal axes depend on the 
first two moments of the distribution. For different ratios 
of the diffusion coefficients we distinguish between oblate 
and prolate cases, which show a concave and a convex 



curvature in the logarithmic representation, respectively. 
Finally, we offer a guide to quantify the eigenvalues of D 
from projected observations and to reconstruct the diffu- 
sion tensor in three dimensions from the moments of the 
tensorial diffusivities. The reconstruction from projected 
observations is possible although any directional infor- 
mation is discarded when determining the distribution 
of diffusivities. 

In summary, the distribution of diffusivities provides 
an advanced analysis of anisotropic diffusion processes. 
The distribution is easily obtained from measured trajec- 
tories and allows for a characterization of the processes. 
Based on the distribution, anisotropic systems can be 
distinguished qualitatively from isotropic processes. Fur- 
thermore, from the moments of the distribution the dif- 
fusion coefficients can be reconstructed. Beyond that, 
the concept of diffusivities as scaled displacements is 
extended to tensorial diffiisivities, which allow the re- 
construction of the diffusion tensor from their first mo- 
ments. Hence, the distribution of diffusivities comple- 
ments well-established methods, such as mean squared 
displacements, for the analysis of diffusion data. 

In future research, we will address the distinction be- 
tween anisotropic and heterogeneous diffusion processes, 
which also involve more than one diffusion coefficient. 
Furthermore, since the eigenvalues of tensor D are in- 
variant to orthogonal transformations, we will apply our 
distribution of diffusivities to systems where the diffusion 
tensor changes its orientation in space and time, such as 
diffusion of ellipsoidal particles in isotropic media and 
diffusion in liquid crystalline systems with an inhomoge- 
neous director field. 
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